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Genome-wide Study of Gene Copy Numbers, 
Transcripts, and Protein Levels in Pairs of 
Non-invasive and Invasive Human Transitional 
Cell Carcinomas* 
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Gain and loss of chromosomal material is characteristic 
of bladder cancer, as well as malignant transformation In 
fleneraL Tho consequences of these changes at both the 
transcription and translation levels Is at present unknown 
partly because of technical limitations. Hem we have at* 
tempted to address this xjuestion in pairs of non-Invasive 
and Invasive human bladder tumors using a combination 
of technology that Included comparative genomic hybrid- 
zatfon, high density oligonucleotide array-based monitor- 
ing of transcript levels (5600 genes), and high resolution 



phenomenon at both the transcription and translation .levels. 
High throughput array studies of the breast cancer cell line 
BT474 has suggested that there Is* a correlation between 
DNAcopy numbers and gene expression in- highly amplified 
areas (2), and studies of Individual genes In solid tumore 
have revealed a good correlation between gene dose and 
mRNA or protein levels In the case of c-erb-B2. cyclln d1 t ' 
ems1 f and N-myc (3-5). However, a high cycfln Dt protein 
expression has been observed without simultaneous am- 
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™^ men fonal gel electrophoreWfhe results showed^^tlon <4), and a low level of c-myc copy number in 
that there Is a gene dosage effect feat In some cases^ crease was observed wirhm* r^^x^ Z^ — 



that there Is a gene dosage effect W In some cases 
superimposes on other regulatory mechanisms. This ef- 
fect depended (p < 0.015) on the magnitude of the com- 
parative genomic hybridization change. In general (18 of 
# ^f^' c ^ ro ^somaf areas with more than 2-fo!d gain 
QT DNA showed a corresponding Increase In mRNA tran- 
scripts. Areas with loss of DNA, on the other hand, 
showed either reduced or unaltered transcript levels} Be- 
cause most proteins resolved by two-dimensional gels 
•are unknown it was only possible to compare mRNA and 
protein alterations In relatively few cases of well focused 
abundant proteins- ?flfrth few exceptions we found a good 
correlation <p < 0.00?) between transcript alterations and 
protein levels. The Implications, as well as limitations, 
of the, approach are discussed. Molecular & Cellular 
PioieomlGs 1&7-45, 200Z 



. Aneuploldy is a common feature of most human cancers 
(1), but little is known about the genome-wide effect of this 
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crease was observed without concomitant c-myc protein 
overexpresslon (6). . 

In human bladder tumors, karyotyping, fluorescent In sftu 
hybridization, and comparative genomic hybrid^zat^on•(CGH) , 
have revealed chromosomal aberrations that seem to be 
characteristic of certain stages of disease progression. In the 
case of non-Invasive pTa transitional cell carcinomas (TCCs), 
this Includes loss of chromosome 9 or parts of it, as well as 
loss of Y In males. In minimally Invasive pT1 TCCs, the fol- 
lowing alterations have been reported: 2q- f 11p-, iq+ . 
1lq13+, 17q+> and 20q+ (7-12). ft has been suggested that 
these regions harbor tumor suppressor genes and onco- 
genes; however, the large chromosomal areas involved often 
contain many genes, making meaningful predictions of the . 
functional consequences of losses and gains very difficult. 

In this Investigation we have combined ge&ome-wide tech- 
nology for detecting genomic gains and losses (CGH) with 
gene expression profiling techniques (mlcroarrays and pro- 
teomlcs) to determine the effect of gene copy number on 
transcript and protein levels In pairs of non-Invasive and Irv 
vaslve human bladder TCCs, 

EXPERIMENTAL PROCEDURES 
MaferfaZ-Bladder tumor biopsies, were sampled after Informed 
consent was obtained and after removal of tissue for routine pathol- 
ogy examination. By light microscopy tumors 335 and 532 were 
staged by an experienced pathologist as pta (superficial paplfcry) 
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The abbreviations used are: CGH, comparatfva genomic hybrid- 
fcatlon; TCG. transitional cell carcinoma: UOH, loss of heterozygositv 
PA-FABP. psoriasis-associated fatty acid-binding protein- 2D 
two-dimensional. ' 
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compared with the n^nv^T^teSa? 2S ^B^^J^SSTT^ expressl ? n of ««*A *» ""ash* tumor 733 as 
■ coua.erpart tamer Sa^TheaveragXSrtSSo ta^KSEr !2™ "T 827 , C ° mpaKKf ^ me Evasive 
• Ikft): The bo/d omih the rattojW raS^n^ md ^ D ^b shown along 0^ length of the c^romoaome 
deviation. The cenfra/ vartfcSZfc^S^ VSS ^u«^m ?^T n ^^r n ^ ed * ,Nn ct " ves tat « caan 8 <™ standard 
OS (/eft) and 2.0 (r&tt). lnXm^^SfrX^!l^l 1 fe,*"^ 80(1 vertka, "n^ next to It (dotted)" Indicate a raBo of 
profile of that dSSZto tt2^«f^SI ^or335 used for comparfeon showed altemtiona In DNA content, tte ralfo 

the geile, and IheloZwSK e^K^o^L fs DLDK/sdataJ,tm1). The bars Indicate the purported tocaft* , of 

Increase (Mae*), >^SSS™,°i ftB ^ * th 1 9 JT vas,ve ,umor «™P^ with the noninvasive counterpart; >2-foW 
InexpresSnXcl^ 

<«erm,a,,^ 



grade I and 0, respectively, tumors 733 end 827 were staged as pTt 
Onve^ve Into wbmucosa). 733 was staged as solid, and 827 was 
staged as papBlary. both grade III. 

J!! fl ^^ a ^° f, - T,65Ue biopsies, obtained fresh from surgery, 
embedded Immediately In a sodium-guanldlnlum mlocyanate 

teolaUo " method (WAK-Cfiemle Medical GMBH). 
Sa IdftSg^ ^ ty " 0li90(dT) M,ectlon **> COUgotex 
•n,° fV f paraff0 '!- 1 W of mRNA was used as starting material 
V^^Tt^T^ e ^ d cONA synthesis was performed using the 
U ^^!l P .f ° ,Ca system Onvttrogen) according to the manufac- 
turer's Instructions but using an oligo(dT) primer containing a T7 RNA 

fflSf? ^2*? SBa Ube,Bd CRNA * WW* «*»9 the ME- 
QAscrlp® /„ w/m transcrtpilon kit (Amblon). BiotfrHabeled CTP and 



OTP (Enzo) was used, together with unlabeled NTPs In the reaction. 
Fdtewlng the In vffro treneedption reaction, the unincorporated nu- 
cleotides were removed using RNeasy columns (Qfagen). 

Amy Hybridization and Scano/ngr-Array hybrtdlzaUon and scam 
ntng was modified from a previous method (13). 10 ng of cRNA was 
fragmented at 94 'C for 3S mln In buffer containing 40 mu Trte 
acetote. pH 8.1, 100 tm KOAc. 30 raM MgOAo. Prior to hybridization 
toe fragmented cRNAlnafix SSPE-T hybridization buffer (1 MNaCf 
10 mu Trfs. pH 7.6, 0.005% Triton), was heated to 85 «C for 5 mm 
subsequently cooled to 40 "C. and loaded onto the Asymetrix probe 
array cartridge. The probe array was then Incubated for 16 h at 40 °C 
at constant rotation (SO rpm). The. probe array was exposed to 10 
washes In 6x SSPE-T at 25 'C followed by 4 washes In O&x SSPE-T 
at 60«C. The btotinylated cRNA was stained, with a streptavldfi,- 
phycoeiythnn conjugate, 10 /xg/ml (Molecular Probes) In 6x SSPE-T 
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Fig. 1— continued 



forSO mln at 25 °C followed by 10 washes la 6x SSPE-T at25 D C.The 
probe arrays were "scanned at 560 nm using a corrfocaJ laser scanning 
microscope (made for Affymetnx by Hewlett-Packard). The readings 
from the quantitative scanning were analyzed by AflyrnetrU gene 
expression analysis software. 

^Mlcn>$steiStQ Ana^^lAiGn^m\\B Analysts was performed as 
described previously (14). MicrosateltUes were selected by use of 
wwwJKJbLr0m.nIh.g6v/genernapg8, and primer sequences were ob- 
tained from the genome data base atwww.gdb.org. DNAwas extracted 
from tumor and blood and amplified by PCRIna volume of 20 pi for 35 
• cycles. The ampllcons were denatured and electrophoresed for 3 h In an 
ABi Prism 377. Data were collected fn the Gene Scan program for 
fragment analysis. Loss of heterozygosity was defined as "fess than 33% 
Of one allele detected In tumor ampflcons compared with Wood. 

Profeomfe Anafyste—TGCa were minced -Into email pieces and 
homogenized In a small glass homogenlzer In 0.5 ml of lysis solution. 
Samples were stored at -20 °C untO use. The procedure for 2D gel 
electrophoresis has been described in detail elsewhere (15, 16). Gels 
.were stained with sflver nitrate and/or Cocrnassle Brilliant Blue. Pro- 
teins were Identified by a combination of procedures that Included 
rn?crosequencIng t mass spectromeuV, two-dimensional gel Western 
ImmunoWottlng, and comparison with/the master two-dimensional gel 
Image of human Keratlnocyte proteins; see Wobase.dk/cgl-Wn/ceIIs, 

CGH— Hybridization of differentially labeled tumor and norma! DNA 
to normal metaphase chromosomes was performed as described 
previously (10), Fluoresceln-labeled tumor DNA (200 ng) f Texas Red- 



labeled reference DNA (200 ng), and human Cot-1 DNA {20 fig) were 
denatured at 37 °C for 5 mln and applied to denatured normal rnet- 
aphase slides. Hybridization was at 37 °C for 2 days. After washing, 
trie sBtfes were oounterstained with 0.15 jtg/ml 4,(^amtdTno-2^phe- 
nyiindote fn an anU-Tade solution. A second hybridization was per- 
formed for aJl tumor samples using fluorosceJn-labeied reference DNA~ 
and Texas Red-labeled tumor DNA (inverse labeling) to' confirm the 
aberrations detected during the initial hybrfcfizatfoa Each CGH ex- 
perfmerft also Included a normal control hybridization using fluores- 
cein- and Texas Red 4ab sled normal DNA. Digital Image analysis was 
used to Identify chromosomal^ regions with abnormal ftjorascence 
ratios/ Indicating; regions of DNA gains and losses. The average 
gmenred fluorescence Intensity ratio profiles were calculated using 
four Images of each chromosome (eight chromosomes total) with 
normaDzaUon of the greervred fluorescence Intensity ratio for the 
entire metaphase and background correction Chromosome Identifi- 
cation was performed based on 4 v 6Hilamld^o-2>phenyllndol8 band- 
ing patterns. Only Images showing uniform high Intensity fluores- 
cence with minimal background staining were analyzed All ' 
centromeres,- p arms of acrocentrio chromosomes, and heterochro- 
matic regions were excluded from the analysis. 

RESULTS 

Comparative Genomic Hybridization- The CGH analysis 
Identified a number, of chromosomal gains" and losses In the 
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Table I 

Correlation between alterations detected by CGH end by expression monitoring 

• Top, CGH used as Independent variable (ff CGH alteration - what expression ratio was found); bottom, altered expression used as 
independent variable (ft expression alteration - what CGH deviation was found). 



CGH alterations 



1\<mor733 vs>335 
Expression change clusters 



Concordance CGH alterations 



Tumor 827 vs. 532 
Expression change clusters 



Concordance 



13 Gain ■ 10 Up-regutatJon 
0 Down-regulation 
$ No change 

• 10 Loss 1 Up-regulation 

6 Down-regulation 
4 No chango 



l 77% 
60% 



tOGaln. 
12 Loss 



6 Up-regulation 
OCovwwegulatlon 

2 No change 

3 Up-regulatlon 

2 Down regulation 

7 No_chango. 



Expression change clusters 
16 Up-regulation 

21 Down-regulation 

15 No change 



Tumor 733 vs. 335 
CGH alterations 



Concordance Expression change dusters 



Tumor 827 vs.532 
CGH alterations 



8096 
1796 

Concordance 



11 Gain 

2 Loss 

3 No change 
1 Gain 

8 Loss 

12 No change 
3 Gain 

3 Loss 

9 No change 



38% 
60% 



17 Up-regulation 
9*DoVYn-regulatlon 
21 No change 



10 Gain 
5 Loss 

2 No change 
OGaln 

3 Loss 

8 No change 
1 Gain 
3 Loss 

17 No change 



5996 
33% 
81 % 



; two Invasive tumors (stage pT1. TCCs 733 and 827), whereas 
the two non-Invasive papillomas {stage pTa, TCCs 335 and 
532) showed onfy 9p-, 9q22-q33- r and and 7+, 9q-, 
and Y-, respectively. Both Invasive tumors showed changes 
(1q22-24+, 2q14.1-qter- F 3q12-q13.3~, 6q12-q22~, 
9q34+ , 11q12-q13+, 17+ T and 20q11.2-q12+) that are typ- 

' leal for their disease stage, as well as.addftional alterations, 
eome of which are shown In Fig. 1. Areas with gains and 
losses deviated from the normal copy humber to some extent, 
and the average numerical deviation from normal was 0.4-fotd 
In the case of TCC 733 and 0.3~fo!d for TOO 827. The largest 
changes, amounting to at least a doubling of chromosomal 
content, were observed at 1q23 In TCC 733 {Fig, 1/4) and 
20q12 In TCC 827 (Rg. 18). 

1 mfiNA Expression In Relation to DNA Copy Number—The 
mRNA levels from the : two Invasive tumors^TCCs 827 and 
733) were compared with the two noninvasive counterparts 
..(TCCs 532 and 335). This was done In two separate experi- 
ments In which we compared TCCs 733 to 335 and 827 to 
532, respectively, using two different scaling settings for the 
arrays to rule but scaling as a confounding parameter, Ap- 
proximately 1,800 genes that yielded a signal on the arrays 
were searched In the Unlgene and Genemap data bases for 
chromosomal location, and those with a known location 
(1096) were plotted as bars covering their purported locus. In 
that way It was possible to construct a graphic presentation of 
DNA copy number and relative mRNA levels along the indi- 
vidual chromosomes (Fig, 1). 

For each mRNA a ratio- was calculated between the level In 
the Invasive versus the non-Invasive counterpart Bars, which 
represent chromosomal location of a gene, were color-coded 
according to the expression ratio, and only differences larger 



than 2-fold were regarded as informative (Fig. 1). The density 
of genes along the chromosomes varied, and areas contain- 
ing only one gene were excluded from the calculations. The 
resolution of the QGH method is very low, and some of the 
outlier data may be because of the fact that the boundaries of 
trie chromosomal aberrations are not known at high resolution. 

Two sets of calculations were made from the data Fbr the 
first set we used CGH alterations as the Independent variable 
and estimated the frequency of expression alterations In these 
chromosomal areas. In general, areas with a strong gain of 
chromosomal material contained a cluster of genes having 
Increased mRNA expression. For example, both chromo- 
somes 1q21-q25, 2p arid 9q, showed e relative gain of more 
than 100% In DNA copy number that was accompanied by 
Increased mRNA expression levels In the two tumor palrs'pg. 
. 1). In most cases, chromosomal gains detected by CGH were 
accompanied by an Increased level of transcripts In both 
TCCs 733 (TT%) and 827 (803*) (Table I, top). Chromosomal 
losses, on the other hand, were not accompanied by de- 
creased expression In several cases, and were often regis- 
tered as having unaltered RNA levels (Table I, top). The IrtablH 
fry to detect RNA expression changes in these cases was not 
because of fewer genes mapping to the lost regions (data not 
shown). 

In the second set of calculations we selected expression 
alterations above 2-fold as the Independent variable and es- 
timated the frequency of CGH alterations In these areas. As. 
above, we found that Increased transcript expression corre- 
lated with gain of chromosomal material (TCC 733, 69% and 
TCC 827, 59%), whereas reduced expression was often de- 
tected In areas with unaltered CGH ratios (Table I, bottom): 
Furthermore, as a control we looked at areas with no after- 
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Expression changes 



Expression changes' 
not detected 



Expression changes 
detected 



Expression changes 
* net delected 



Tumor 827 versus 532 Tumor 733 versus 335 

•n^r! 0 ^ 3 ! 0 " b j il T* maxImurn CGH aberration and the ability to detect expression change by oligonucleotide emiy 
SSS^^W Z!!* 0 ™ 1 95 *. numerfcaI * " change In ratio between Invasive tumors 827 (A) and 733 (♦ ) arid their noninvasive 
^^fj? 1 ^expression change was taken from the Expression Gne to the right In Hg. 1, which depicts the resulting 
expr^pnehange for a given chromosomal region. At feast half of the mRNAs from a given region have to be either up- or dovwweguiated 
rat? value*!* 5 ^ T" 9 ** M Chmm050mal "«* 1(1 VJbkh TOH mUo * Ius * m,nl * °™ standard deviation was ouSdTthe' 



atlon In expression. No alteration was detected, by CGH In. 
most of these areas (TCC 733, 60% and TCC 827, 81 %;see 
Table I, bottom). Because the ability to observe reduced or 
Increased mRNA expression clustering" to a certain chromo- 
.sornal area clearly reflected the extent of copy number 
changes, we plotted the maximum CGH aberrations In the 
regions showing CGH changes against the ability to detect a 
change In mRNA expression as monitored by the oligonucleo- 
tide arrays (Rg. 2K£or both tumors TCC 733 fp < 0.015] and 
TCC 827 (o < 0.00003) a highly significant correlation was 
observed between the level of CGH ratio charge (reflecting 
the DMA copy number) and alterations detected by the array 
based technology (Rg. 2^ Similar data were obtained when 
areas with altered expression were used as Independent vari- 
ables. These areas correlated best with CGH when the CGH 
ratio deviated 1.^ to 2.0-fold (Table I, bottom) but mostly did 
not aj lower CGH deviations. These data probably reflect that 
loss of an allele may only lead to a 50% reduction In expres- 
sion level, which Is at the cut-off point for detection of expres- 
sion alterations. Gain of chromosomal material can occur to* a 
much larger extent 

Mlctosatefflte-based Detection of Minor Areas of Loss- 
es—In TCC 733, several chromosomal areas exhibiting DNA 
amplification were preceded or followed by areas with a nor- 
mal CGH but reduced mRNA expression (see Rg. 1 , TCC 733 
chromosome 1q32, 2p21, and 7q21 and q32, 9q34, and 
10q22). To determlne whether these results were because of 
undetected loss of chromosomal material in these regions or 



because of other non-structural mechanisms regulating tran- 
scription, we examined two mlcrosatellltes positioned at chro- 
mosome 1q26-32 and two at chromosome 2p22. Loss of 
heterozygosity (LOH) was found at both 1q25 and at 2p22 
Indicating that minor deleted areas were not detected with the 
resolution of CGH (Rg. 3). Additionally, chromosome 2p In 
TCC 733 showed a CGH pattern of gain/ho change/gain of 
DNA that correlated with transcript Increase/decrease/Ir^ 
crease. Thus, for the areas showing Increased expression 
there was a correlation with the DNA copy number alterations 
(Rg' 1/4). As Indicated above, the mRNA decrease observed In 
the middle of the chromosomal gain was because of LOH, 
Implying that one of the mecHanlsms for mRNA down-regu- 
lation may be regions that have undergone smaller losses of 
' chromosomal material. However, this cannot be detected with 
the resolution of the CGH method;; " " — 

In both TCC 733 and TCC 827, the telomertc end of chro- 
mosome 11 p showed a normal ratio In the CGH analysis; 
however, clusters of five and three genes, respectively, lost 
their expression. Two mlcrosatellltes (D11S1760, 011S922) 
positioned dose to MUC2, IGF2, and cathepsln D Indicated 
LOH as the most likely mechanism behind the loss.of expres- 
sion (data not shown). 

A reduced expression of mRNA observed in TCC 733 at 
chromosomes 3q24 t 11p11, 12p12.2. 12q21.1, and 16q24 
and in TCC 827 at chromosome 11p15.5 f 12p11, 15q11.2, 
and 18q12 was also examined for chromosomal losses using " 
mlcrosatellltes positioned as close as possible to the gene loci 
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^Fkx 5. Miwosataliite analysis of loss of heterozygosity. Tumor 
733sho5jrtng toss of heterozygosity at chromosome 1q25, detected 
W by D1S21^ctose to Hu class I hlstbcompaUbllHy antigen (gene 

. vm ^^ *fl* 1 >> W !>y D1S2735 dose to cathepsln E (gene 
numbers In Rg. 1), and (c) at chlomosome 2p23 by D2S2251 close 

• to general ^-spectrin (gene number 1 1 on Rg. 1) and of (d) tumor 627 
snowing toss of heterozygosity at chromosome 18q12 by S18S1118 
<%>f*^ nrilrochofKlriaf 3-oxoacyI~coenzyme A thiolase (gene numper 
«?• 1)? The upper cv/ves show the efectropherogram obtained 
from normal DNA from leukocytes (/V). and the tower curves show the 
efectropherogram from tumor DNA (7). In all cases one allele Is 
partteJfy lost fan the tumor ampllcon. 

showing reduced mRNA transcripts. Only the mlcrosatelHte 
positioned at 18ql2 showed LOH (Fig. 3), suggesting that 
tmnscrifitlonaJ.down-regula«on of genes In the other regions 
may be controlled by other mechanisms. 

Relation between Changes In mRNA and Protein Levels- 
2D-PAGE analysis. In combination with Coomassle Brilliant 
Blue and/or silver staining, was canied out on ail four tumors 
using fresh biopsy material, 40 well resolved abundant known 
proteins migrating In areas away from the edges of the pH 



E 1 



It 



"A" 



B 

4W- 



fteduced 
protein 



Increased 
protein 



Ra 4. Correlation between protein levels as Judged by 20- 
PAGE and transcript ratio. For comparison proteins were divided In 
Ihree groups, unaltered In level or up-or down-regulated {horizontal 
axis). The mRNA ratio as determined by oligonucleotide arrays was 
plotted tor each gene {vertical wets), a,; mRNAsjhat were scored as 
present In both tumors used for the ratio calculation; A, mRNAs that 
were scored as absent In the Invasive tumors {along rwrfeonfa/a»s) or 
as absent In noninvasive reference {fop of figure). Two different 
scallngs were used to exclude scaling as a confounder; TCCs 827 
and 532 (k&) were scajed with background suppression, and TCCs 
733 and 335 were scaled without suppression." Both compari- 
sons showed highly significant {p < 0.005) differenceeln mRNA ratios 
between the groups. Proteins shown were as follows: Group A {from 
/efl), phosphoglucomutftse 1, glutathione transferase class & number 
4, fatty acid-binding protein homologue, cytokeratln 15, and cyto- 
keratin 13: B {from fefl), fatty acld-blndlng protein homologue, 28-KDa 
heat shock protein, cytokeratln t3 f and calcyclln; C<firom tefl), *-enc~ 
lase, hnRNP B1, 28-kDa heat shock protein, 14-3-3-e, and 
pre-mRNA splicing factor; O, mesoiheJlBl keratin K7 (type H); £ (from 
top), glutathione S-transferase-Tr and mesotheflal keratin K7 (ly^o.ll)' 
F(from top and feft), adenyfyt cyetase-assoclated protein, Bcadherin 
keratin 19, calgtaartn, phosphogtycerate mutasc.annexln IV, cy-- 
toeketeta! *actfn. _ho.RNR. A1,. Integral membrane protein caJnextn 
(IP90), hnRNP H t brain-type dathrln light chaln-a, hnRNP F t 70-kDa 
heat shock protein, heterogeneous nuclear rfronudeoprotetn A/B 
translatlonalty controlled tumor protein, liver gtyceraldehyde-O-phos^ 
phate dehydrogenase, keratin 8, aldehyde .reductase, and Na,K- - 
ATPase p-1 subuntt; Q. (from to# and feff), TCP2Q, calgfezarla 70- 
kDa heat shock protein, calnexfn, hnRNP H, cytokeratln 15, ATP 
synthase, keratin 1 9, .triosephosphate Isome rase, hnRNP F. Rver gjyc- 
eraIdehyde-3-phosphatase dehydrogenase, glutathione S-UBnsfer-7 
ase-ir, and keratin 8; H (from feft), plasma gBfedfln, autoantlgen cal- 
reticuRn, thiorodoxJn, and NAD+-dependent IShyo^vpn^iajidln 
. dehydrogenase; f (from top), prolyl 4-hydroxylase 0-eubunH, cyto- 
keratin 20, cytokeratln 17, prohibition, and fructose ^^phos- 
phatase; J annexln II; K, annaxln IV; L (from fop and fefl), 90-kDa heat, 
ehock protein, prolyl hydroxylase p-eubuntt, o-enolase, GRp 78, 
cyclophllrn, and cofOln. 

gradient, and having a known chromosomal location, were 
selected for analysis In the TCC pair 827/532. Proteins were 
Identified by a combination of methods (see "Experimental 
Procedures 0 ). In general there was a highly significant corre- 
lation (p < 0.005) between mRNA and protein alterations (Rg. 
4). Only one gene showed disagreement between transcript ~ 
alteration and protein alteration. Except for a group ofcyto- 
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Fie. 5. Comparison of protein and transcript levels In Invasive 
andnon-IriVBstve TCCs. The upper part of the figure shows a 2D gel 
(fefl) and the oligonucleotide array (rfcM) of TCC 632. "The red lecfan- 
glas on the upper ge/ highlight the areas that are compared below. 
Identical areas of 2D gels of TCCs 532 and 827 are shown below. 
Clearly, cytokerattns 13 and 15 are strongry down-regulated In TCC 
827 (ped annotation). The tne on the array containing probes for 
cytokemtfn.15 Is enlarged below the array (reef arrow) from TCC 632 
and Is compared with TCC 827. The upper row of squares In each tile 
corresponds to perfect match probes; the tower row corresponds to 
mismatch probes containing a mutation (used for correction for un- 
speeffic binding). Absence of signal Is depicted as black, and the 
higher the signal the fighter the color, A high transcript level was 
detected In TCC 532 (6151 units) whereas a much lower level was 
.detected In TCC 827 (absence of signals), for cytokeratln 13, a high 
transcript level was also present In TCC 532 (15659 units), and a 
much lower level was present In TCC 827 (623 units). The 2D gels at 
the Oottom of the flgure (fen} show levels of .PA-FABP and adipocyte- 
FABP In TCCs 336 and 733 (Invasive), respectively. Both proteins are 
down-regulated In the Invasive tumor. To the right we Show the array 
tiles for the PA-FABP transcript. A medium transcript level was de- 
tected In ths-case of TCC 335 (1*277 unite) whereas very low levels 
were detected In TCC 733 (166 units). l&. Isoelectric focusing. 



keratins encoded by genes on chromosome 17(Rg, 5) -the 
analyzed proteins did not belong to a pajtteular'famliy* 2S well 
focused proteins whose genes had a know chromosomal 
location were detected In TCCs 733 and 335, and of these 1 9 
correlated {p < 0.006) with the mRNA changes detected using 
the arrays (Fig, 4). For example, PA-FABP was highly ex- 
pressed In the non-Invasive TCC 335 but lost In the Invasive 
counterpart (TCC 733; see Fig, 5). The smaller number of 
proteins detected In both 733 and 335 Was because of the * 
smaller sfze of the biopsies that were available. 

11 chromosomal regions where CGH showed aberrations 
that corresponded to the changes In transcript levels also 
showed corresponding changes In the protein level (Table II). 
These regions Included genes that encode proteins that are 
found to &e frequently altered In bladder cancer, namely 
cytokeratins 17 and 20, annextns I! and IV, and the fatty 
acld-blndlng proteins PA-FABP and FBPl! Fbar of these pro- 
teins were encoded by genes In .chromosome T7q t a fre- 
quently amplified chromosomal area In Invasive bladder 
cancers. 

DISCUSSION 

Most human cancers have abnormal DNA content, having 
lost some chromosomal parts and gained others. The present 
study provides some evidence as to the effect of these gains 
and losses on gene expression In two pairs of non-lnvasjve 
and Invasive TCCs using high throughput expression arrays 
and proteomics, In combination with CGH. In general, the 
results showed that there Is a clear Individual regulation of the 
mRNA expression of single genes, which In some cases was 
superimposed by a DNA copy number effect Inmost cases, 
genes located In chromosomal areas with gains often exhib- 
ited Increased mRNA expression, whereas areas showing 
losses showed either no change or a reduced mRNA expres- 
sion. The latter rrjlght be because of the fact that losses most 
often are restricted to loss of one allele, and the cut-off point 
for detection of expression alterations-was a 2-fbkf change.- 
thus being at the border of.detection. In several cases, how- 



\Table II 

Proteins whose expression level correlate? with both mRNA and gene dose changes 



Protein 



Chromosomal location Tumor TCC CGH alteration Transcript alteration" Protein alteration 



Annexf n II 
Annexln IY 
Cytokeratln 17 
Cytokeratln 20 
(PA-JFABP . . 

FBPT 

Plasma getsolin 
Heat shock protein 28 
Prohibffln 
ProtyM-hydroxyl 
hnRNPBI 



1q21 


733 


Gain 


AbstoPres 0 


Increase 


2p13 


733 


Gain 


3,9-Fbtd up 


Increase 


' 17q12-q21 


827 


Gain - 


3.8-Fold up 


Increase 


17q21>1 


827 


Gain 


5,6-Fold up 


Increase 


SqZI.2 


827 


loss 


10-Fold down 


Decrease 


9q22 


827 


Gain 


2,3-Fold up 


Increase 


9q31 


827 


Gain 


Absto Pres 


■Increase 


15q12-q13 


827 


Loss 


2.5-Fold up 


Decrease 


17q21 


827/733 


Gain 


3.7-/2.S-Fold up 6 


Increase 


17q25 


827/733 


Gain ' 


6.7V1,e-Fotdup 


increase 


7p15 


827 


Loss 


2.5-FoW down 


. Decrease 



°Abs, absent; Pres, present 

* In cases where Ihe correspondfng alterations were found in both TCCs 827 and 733 these are shown as 827/733. 
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ever, an Increase or decrease In DMA copy number was 
associated with cte novo occurrence or complete loss'.of tran- 
script, respectively. Some of these transcripts could not be 

. detected In the non-Invasive tumor but were present at rela- 
tively high levels In areas, with DNA amplifications In the Inva- 
slveturrwrsfeg. InTCC 733 transcript from'cellular Ilgand of 
annexion gene (chromosome 1q21) from absent to 2670* 
arbitrary Units; In TCC 827 transcript from small proline-rich 
protein 1 gene (chromosome 1q12-q21.1) from absent to 
1326 arbitrary units). It may be anticipated from these data 

_that significant clustering of "genes with an Increased expres- 
sion to a certain chromosomal area Indicates an Increased 
likelihood of gain of chromosomal material In this area. 

Considering the many possible regulatory mechanises act- 
ing at the level of transcription, it seems striking that the gene 
dose effects were so clearly detectable In gained areas. One 
hypothetical explanation may lie In the loss of controlled 
methylation In tumor cells (17-19). Thus, It may be possible . 

•that la chromosomes with Increased DNA copy numbers two 

. or more alleles could be demethyiated simultaneously leading 
to a higher transcription level, whereas In chromosomes with 
tosses the remaining allele could be partly methylated, turning 
off the process (20, 21). A recent report has documented a 
ptoldy regulation of gene expression In yeast, but In this case all 
the genes were present In the same ratio (22), a situation that is 
not analogous to that of cancer cells, which show marked 
chromosomal aberrations, as well as gene' dosage effects,- 

Several CGH studies of bladder cancer have shown that 
some chromosomal aberrations are common at certain 
stages of disease progression, often occurring in more than 1 
of 3 tumora In pTa tumors,.these Include 9p-, 9q-, 1q+> Y— 
(2, 6), and ln4>T1 tumors, 2q-,11p- l 11q~, 1q+, 5p+, 8q+, 
17q+, and 20q+ (2-4, 6, 7). The pTa tumors studied here 
showed similar aberrations such as 9p- and 9q22-q33- and 
9q— and respectively, Ukewfee, the two minimal invasive 
pT1 tumors showed .aberrations that are commonly seen at 

Jfaatitoge^andTC^^ 
commonly seen pattern of losses and gains, such as 1q22-24 

* amplification (seen In both tumors), 1 1q14-q22 loss, the latter 

. often linked to 17 q+ (both tumors), and 1q+ and 9p~, often 
. linked to20q+ and 11 q13+ (both tumors) (7-9). These ob- 
setvatlons.lndicate that the pairs of tumors used In this study 
exhibit chromosomal changes observed In many tumors, and 
therefore the findings could be of general. Importance for 
bladder cancer. 

Considering that the mapping resolution of CGH Is of about 
•20 megabases It is only possible to get a crude picture of 
chromosomal Instability using this technique. Occasionally, 
we observed reduced transcript levels dose to or Inside re- 
gions with Increased copy numbers. Analysis of these regions 
by positioning heterozygous microsatellltes as close as pos- 
sible to the locus showing reduced gene expression revealed 

' Ic&s of heterozygosity In several cases. It-seems likely that 
multiple and different events occur along each chromosomal 



arm and that the use ofcDNA mlcroarraysfor analysis of DNA 
copy number changes will reach a resolution that carl resolve 
these changes, as has recently been proposed (2). The outlier 
• data were not more frequent at the boundaries of the CGH 
aberrations. At present we do not know the mechanism be- 
hind .chromosomal aneuplotdy and cannot predict whether 
chromosomal gains will be transcribed to a larger extent than 
the two native alleles. A mechanism as genetic imprinting has 
an Impact on the expression (eve! In normal cells and Is often 
reduced In tumors. However, the relation between Imprlming 

and g^ln of chromosomal material Is not known* * • . • 

We.regard It as a strength of this Investigation that we were 
able to compare invasive tumors to benign tumors rather than 
to nor/nai urothellum, as the tumors studied were biologically 
very close, and probably may represent successive steps in 
the progression of bladder cancer. Despite the limited amount 
of fresh tissue avail able It was possible to apply three different 
state of the art methods. The observed correlation between 
DNA copy number and mRNA expression Is remarkable, when 
one considers that different pieces of the tumor biopsies were 
used for the different sets of experiments. This Indicate that 
bladder tumors are relatively homogenous, a notion recently 
supported by CGH and LOH data that showed a remarkable 
similarity even between tumors and distant metastasis (1 0, 23). 

.In the few cases analyzed, .mRNA. and protein levels 
showed a striking correspondence although in some cases 
we found discrepancies that may be attributed to translations 
regulation, post-translational processing, protein degrada- 
tion, or a combination of these. Some transcripts belong to 
undertranstated mRNA pools, which are associated with few 
translationalry Inactive rlbosomes; these pools, however, 
seem to be rare (24). Protein degradation, for example, may 
be very important. in the case of polypeptides with a short 
half-life (e.o;.-sIgnaling proteins). A poor correlation between 
mRNA and protein levels was found in liver celts as deter- 
mined by arrays and 20-PAGE (25), and a moderate correla- 
tion wasrecentty reported by Ideker etalJgS) in yeast 



Qnterestingiy, cur study revealed a much better correlation 
.between gained chromosomal areas and Increased mRNA 
levels than between loss of chromosomal areas and reduced 
mRNA levels. In general, the level of CGH change determined 
the ability to 'detect a change In transcript) One possible 
explanation could be that by losing one allele the change In 
mRNA level Is not so dramatid as compared with gain of 
materia), which can be rather unlimited arid may lead to a 
severalfold Increase In gene copy number resulting In a much 
higher Impact on transcript level. The latter would be much 
easier to detect on the expression arrays as the cut-off point 
was placed at a 2-fold level so as not to be biased by noise on 
the array. Construction of arrays with a better signal to noise 
ratio may in the future .allow detection of lesser than 2-fold 
alterations I n transcript levels, a feature that may facilitate the 
analysts of the effect of loss of chromosomal areas on tran- 
script levels. 
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la eleven cases we found a significant correlation between 
DMA copy number, mRNA expression; and protein level Four 
of these proteins were encoded by genes located at a fre- 
quency amplified area In chromosome 17q.. Whether DNA 
copy number Is one of the mechanisms behind alteration of 
these eleven protelnsls at present unknown and will have to 
be proved by other methods using a larger number of sam- 
ples. One factor making such studies complicated Is the large • 
extent' of- protein modification that occurs after translation, 
requiring Immunoldentlflcatlon and/or mass spectrometry to 
^coirectiy Mentlfy the proteins In the gels; 

In contusion, the results presented In this study exemplify *, 
the large body of knowledge that may be possible to gather in 
the future by combining state of the art techniques that follow 
the 'pathway from.DNA to protein (26). Here, we used a. tradi- 
tional chromosomal CGH method, but lathe future high reso- 
lution CGH based on mlcroarrays with many thousand radiation 
hybrid-mapped genes' will Increase the resolution and Informa- 
tion dedvedfrom these types of expenments (2). Combined with * 
expression arrays analyzing transcripts derived from genes with 
known locations, and 2D gel analysis to obtain Information at 
the post-translatlortal level, a clearer and more developed un- 
derstanding of the tumor genome will be forthcoming. 
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